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FOREWORD 


This Is an Interim report on the three-dimensional analysis of hulk 
cavitation and secondary water hammer which is being conducted for the Office 
of Naval Research under contract NONR-3709(00). 

In work already carried out by one of the authors for the Office of 
Naval Research, there has been conducted and reported a study of the mechanics 
or physics of bulk cavitation for the one-dimensional case; the work now reported 
is the three-dimensional embellishment. 

We find that, for practicable ranges of yield and depth of burst, the 
size, location, and shape of the cavitated region appears to be quite different 
from that which has been traditionally accepted. 

Thi portion of the analysis bearing on closure and secondary waves is 
still in process, and so is not reported herein; however, as a preliminary, we 
do report on the rudimentary analysis of the trajectory of the water hammer. 

We also make a general discussion of the secondary wave configurations which 
one might expect to be associated with such a traveling hammer. 

We are particularly appreciative of the assistance to our theoretical 
effort which has been provided by Dr. Heinrich Schauer and his associates at 
the Underwater Explosives Research Division—especially the help provided by 
a preliminary analysis of their very recent bulk cavitation experiments. Past 
widespread observations of bulk cavitation phenomena have been most helpful 
throughout this investigation, and we are most grateful for the assistance in 



this regard which has continually been provided by Dr. William Murray of the 
David Taylor Model Basin and Dr. E. H. Swift, Jr., of the Naval Ordnance Labor¬ 
atory. We are also grateful for the support and continual assistance provided 
by the project scientific officer, Mr. James Winchester of the Office of Naval 
Research. 


Respectfully submitted, 
ENC^EERING-PHYSICS COMPANY. 
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THREE-DIMENSIONAL ANALYSIS OF 


BULK CAVITATION — INTERIM REPORT 

I. INTRODUCTION 

This is an interim technical report of work conducted under Contract 
N0NR-3709(00). The work now being conducted is a three-dimensional embellish¬ 
ment of previous work conducted, the final report of which has already been 
submitted. 3, The overall subject matter has consisted of an investigation of 
the underlying theory as well as some of the practicable aspects associated 
with bulk cavitation and subsequent secondary pressure waves which have been 
observed in connection with underwater detonations of conventional high explo¬ 
sives as well as nuclear explosives. 

To date in the three-dimensional effort we have completed, and report 
herein, the detailed analysis—together with machine computational results— 
bearing on the extent of the cavltated region as a function of depth of burst 
and yield. The rigorous or detailed analysis bearing on closure and secondary 
waves is still in process, and so is not discussed in this report; however, as 
a preliminary to the detailed closure analysis, we have carried out and report 
herein a rudimentary analysis of the expected water hammer and its associated 
secondary wave configuration. 

In particular in the work completed so far we have found that the ex¬ 
tent of the cavltated region is quite different in size, location, and shape from 
that which had been traditionally accepted. There has never been an explicit 


Vincent Cushing, "Study of Bulk Cavitation and Consequent Water Hammer," 
Fined Report to Contract NONR-3389(00) for Office of Naval Research, October 31 , 
1961 . 
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* 

bulk cavitation experimental investigation, and so it is not surprising that 
the rather casual indirect observations of bulk cavitation have not portrayed 
either the extent or variety of cavitation which has been found in our models. 

Bulk oavitation is a phenomenon similar to that of spallation; it occurs 
whenever a strong compression wave reflects from a free surface. The reflected 
wave, a rarefaction wave, can give rise to large net tensions throughout a large 
region of the water; since seawater cannot withstand such large tensions the wa¬ 
ter parts or ruptures, spalling layers of water upward. Later, under the influ¬ 
ence of gravitational and other forces, this spalled water again falls back and 
impacts with the underlying unruptured water. This impact or reloading is tanta¬ 
mount to a large water hammer, which sets up strong secondary waves. Indirect 
observation, in connection with past underwater explosion tests, have indicated 
that these secondary waves probably play a large role in the airblast associated 
with such explosions, and also play a large role in the damage sustained by an 
underwater or surface vessel. 

Our analysis, substantiated by reassessment of past field test results, 
indicates that a very large fraction (as high as 70-90 per cent) of the total 
hydrodynamic yield of an underwater explosion is (at least momentarily) absorbed 
in the cavitation process—and this absorbed energy is later re-emitted in the 
water hammer at the time of reloading. This absorption and re-emission process 
transforms completely the character (including wave shape and wave patterns) of 



* 

With the exception of a test conducted in June, 1962 by the Underwater 
Explosives Research Division of the David Taylor Model Basin. At the time of 
writing, data reduction and analysis on this teat data is just commencing, and 
so no results or conclusions are available now. 
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the compression or shock wave originating from an underwater explosion. 

The overall cavltated region does not close everywhere simultaneously. 
Actually, closure and water hammer first occur (under most conditions) along the 
circumference of a circle, the (horizontal) radius of which is a function of yield 
and depth of burst; thereafter the impact point or water hammer spreads horizon¬ 
tally outward and horizontally inward from this so-called ring of first impact. 

In effect, then, we have a traveling water hammer configuration, and this travel¬ 
ing water hammer gives rise to secondary wave patterns. The net effect of this 
traveling water hammer is that the generated secondary wave is generally focused 
toward two points lying on the vertical axis passing through the explosion: one 
focal point appears to lie well beneath the burst point, whereas the second focal 
point lies above the water surface, and thus the reloading secondary wave pattern 
can give rise to large air pressures in the region generally above surface zero. 

In this report we indicate the machine computational results based on a 
rigorous analysis, which describes the overall extent of cavitation as a function 
of depth of burst and yield. In the continuing work, we have set up the analysis 
and are carrying out computations to describe the motion of the cavity boundaries 
during the closure process; and ultimately describe the water hammer and second¬ 
ary wave pattern set up at the time of final closure. 
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II. ANALYTICAL MODEL 


For this study, it is assumed that an explosion takes place underwater 
and that the initial disturbance is a function of radial distance from the 
source (and not a function of angular position). For a small disturbance and 
an adiabatic expansion, the equation of motion is 


5v 1 grad P 

55 " ' 


(i) 


and the equation of continuity is 

1 dP 
“25t 
c o 


■ - p div v 
o 


( 2 ) 


where v » velocity of the fluid at any point; 

p Q = density of the fluid; 

c ■ speed of sound in the fluid; 
o ’ 

P ■ pressure in the fluid at any point. 

For the case to be considered, it is convenient to use spherical coordinates 
with origin at the source, and the assumptions expressed mathematically become 


r dP dP A 

v - u r F ; 55 - 5$ - 0 


(3,4) 


Under these conditions, the equations of motion and continuity become 


* u r 

St" 


_ 1 _ 


5p 

57 


(5) 


1 5p 

~2 St B - p 0 JB 57 


l 5 ,2 v 
( r V 


( 6 ) 
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If u r , the particle velocity, is eliminated, then 


1 5 / 2 dPv 1 d 2 ? 

7 * (r ^ •7 ff 5? 

o 


Any function of the form, 


(7) 


P(r,t) -P -;f (t - r/c ) 


or 

is a solution of this equation. 

If this result is substituted into *uhe equation of motion, then 

t 


u (t) - u (t ) 
r' 7 r' o' 


P - P 


P c 
o o 


[P(r,f) - P 0 ]«. . 


( 8 ) 


(9) 


At this point, it can be observed tl .t, although the fonn of the 

spherical wave does not change as it spreads out, its amplitude decreases due 

to the factor p Also, the particle velocity does not become zero when the 

pressure returns to the equilibrium value, but will continue to have a positive 

value due to the integral term. The development so far has directly followed 
2 

that of Cole. 

Now, if the wave is assumed to be. a shock front with exponential decay, 
f(t - r/c n ) = Ja e” 15 ^ " r / c o } ; t - r/c Q 2 0 
0 ; t - r/c Q < 0 


, ( 10 ) 


‘TRobert H. Cole, Underwater Explosions, Princeton Univ. Press, 
Princeton, N. J., 1948. 
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and 


P(r,t) -P„- ;t-r/c o *0 

J t - r/c < 0 


( 11 ) 


If t a t Is assumed to be a time before the explosion, then u (t ) = 0 and 

0 Y 0 


u r (t) . I jJL. e - b(t - r /=o> ♦ -V [1 - •-*<* • '/«»)) J t - r/c 0 * 0 
*o o o.rT) 


( 12 ) 


0 ; t - r/c Q < 0 

It can be seen that b corresponds to — and — corresponds to — as used 

T C A 

3 ® 

in the one-dimensional analysis already reported. Furthermore, P(r,t) - P q 
corresponds to P 1 in the earlier report, or, in other words, P 1 represented the 
Increase in pressure above equilibrium. Nov, if the depth of the explosion is 
d, then A = P g d, then where P £ is the maximum increase of pressure above equili¬ 
brium at the point on the surface directly above the explosion. If these sub¬ 
stitutions are made, the above equations become 


Pl =PL e - (t / T " r / X > 


; t/r - r/X 2 0 
; t/T - r/X < 0 


, (13) 


“r l-znr «' (t/T ’ rA) + [1 - e- (t/T - r/X) ];t/r - r/X * 0 

° ° P o r . (14) 


;t/ T - r/X < 0 


^Vincent Cushing, loc. 


cit. 
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Now, as in the previous work, let T - t/T, R ■ r/X and P = P/p q c o and, 
in addition, let d/X ■ D, then 



I-RiO 

, (15) 

T - R < O 


u 




T-RiO 

. (16) 

T - R < 0 


One essential difference between the three-dimensional case and the one- 
dimensional case must be noted, T is equal to zero at the time at which the 
shock front would have left the source had it traveled at sonic velocity through¬ 
out its path in this case, whereas, in the one-dimensional case T was equal to 
zero at the time at which the front reached the surface. 

When the disturbance reaches the free surface, it is reflected back as a 
rarefaction. However, in the three-dimensional case, the front does not 
approach the surface normally in general. The reflection is, therefore, not, 
in general, along the same path as the initial wave, but along a line forming 
an angle of the same magnitude with the normal to the surface as the impinging 
disturbance but with opposite sign. The reflection thus follows the familiar 
Snell's Law of Reflection. If the surface is assumed to be a horizontal 
straight line, the reflected wave can be thought of as emanating from a point 
on the same vertical line as the original source at a distance, D, above the 
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surface, the so-called 'toirror image ' 1 of the source. If R- is distance from 
the imaginary or virtual source and a the angle between a given ray and the 
vertical, then Pg, the pressure or tension (since it is a negative pressure) 
in the reflected wave, is given by 



T - R. 4 0, and 

R 1 coS a4 D • (17) 

otherwise 


The geometry of this situation Is sketched in Fig. 1. 
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Fig. 1—Sketch of underwater explosion illustrating 
force field. 
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III. PRESSURE FIELD WITHOUT CAVITATION 


Based on the model developed in Section II above, the total pressure 


at any point within the fluid can be determined as the algebraic sum of 


1 . the compression wave emanating from the source of the explosion; 


2 . the rarefaction wave caused by the reflection of this pressure 


wave at the surface of the fluid; 


3 . hydrostatic pressure; 


4. atmospheric pressure. 


If a coordinate system is set up in the fluid such that the z coordinate it. 


measured downward from the surface and the h coordinate is measured horizon¬ 


tally from the point of the explosion (since it is-assumed that the explosion 


is spherically symmetric, there is no need to have but a single horizontal 


coordinate as the pressure field will be cylindrical^ symnetrlc around the 


vertical line through the explosion), the total pressure at any point (z,h) at 


any time, t, is given by 


gz + P 


t < r/c 


P t (**>,t) 


P s d -(t - r/c )/t . 

e V ' o ,f + gz + P 

* B 


r/c S t < r. /c 
o' o 1' c 




2 2 2 
where r ■ (d - z) + h ; 


P p O 

Tj* - (d + z r + H i 


t is the atmospheric pressure; 
& 
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g is the weight per unit volume of the fluid; and 

P . d, and c are as defined in the previous sect on. 

8 0 

The total pressure may he expressed in a dimensionless form and it is 

in this form that the calculations completed so far have been made. To obtain 

the dimensionless form, each pressure is divided by P , each length is divided 

s 

by t c q and the time is divided by t /. Equation (1) then becomes 


P + (*,h, 



r-(t - r) .-(t - r,) 


e"' 1 - V ‘ 

" ? 1 -I 


+ Gz + P_ 


; t < r 

; r S t < r 3 

; tir x 


(19) 


where each of the quantities is now dimensionless. T^e only tricky constant 
is the G which is equal to the g times r c divided by P . 

O o 

Since this study is concerned primarily with cavitation arid since it 
is known that the presence of cavitation is caused by a negative pressure or 
tension, the time of particular interest at any point is that time at' which 
the rarefaction wave is at its maximum; since the wave being considered has a 
shock front and exponential tailoff, this is the time at which the rarefaction 
front has Just reached the point being considered. This time, of course, in 
the dimensionless presentation of the equation is the time at which t equals r^. 
The pressure at this time at any point (z,h) is given by 

«“ - (t - r) 

P t (z,t) = d[| 1 - i J + Gz + P a , (20) 


where all the qjuscritifcdiari are in dimensionless form. 


W ! 3 * 1 ,vr 
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Equation (20), then, expresses a function which must he studied to 
determine the presence of cavitation although it must he noted that the func¬ 
tion itself. is derived op the tap it assumption that cavitation does not exist. 
Cavitation will begin when P. becomes negative and equ/al to a negative pres- 
sure P^, the breaking or rupture pressure. However, .it must he remembered 
that the rarefaction wave travels along a ray emanating from a "virtual" source 
above the water. .Therefore we. wish to study this function as it varies with 
r^ or along a given ray from the virtual source. 

Representative curves showing the variation of P^ with 
Figs. 2, 3> - and 4. - It can he seen In Fig. 3 that, if the ray being studied 
passes sufficiently near the actual source of the explosion,P t is positive at 
the surface, becomes negative, reaching a minimum, and rises to a large posi¬ 
tive value at or near the explosion, then deoreaseB perhaps becoming negative 
again reaching a mini naan and then rising as r^ increases unt il eventually the 
hydrostatic pressure becomes the dominant term. On the other hand if a ray 
further removed from the explosion is followed, the value of is still slightly 
positive at the surface then becomes negative reaching a minimum, and theh 
begins to ri^e becoming positive and approaching the hydrostatic pressure as r^ 
Increases. It is perhaps obvious that if a ray so far removed that no effects 
of the explosion could be discernible were followed then would begin at its 
surface value and Just increase due to hydrostatic pressure. 

In order to determine horizontal extent for the cavitated region the 
expression was expanded in power series assuming that h was much greater 
than d and z and the point at which the derivative of P^ with respect to r^ 


r^ are 


shljwn 


in 
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was equal to zero- at the surface vas determined through a first order approxi¬ 
mation* , It was found.that this condition is reached when 

h * 3 I 4 O; d 1 / 2 ' ,' (21) 

where h and d are in feet. It is rather interesting to oh serve that the value 
'of h so determined is independent of the yield of the explosion although it 
must be remembered that this is only a first order approximation and is not 
valid for d > 1150 feet, approximately. 

A breakdown of the total pressure into its components appeared to be of 
some interest; so in Fig. 5 the components of one of the curves plotted in' 

Fig. 3 are plotted as a function of r.. P is constant and has such a small 

i - - la. 1 

value (10" J ) that it could not be plotted. It can be seen that the hydrostatic 
pressure is uniformly increasing with r v and the reflected rarefaction pressure 
is uniformly decreasing with r^. However, as we go along a ray from the 
"virtual" source the pressure from the original explosion behaves very errati¬ 
cally and it is the variation of this pressure that produces the unusual facets 
in the total pressure. To look into the source of this variation it was 
thought to be of interest to determine the phase of the compression wave as 
sampled along a given ray at the head of the rarefaction wave. Phase, in 
this sense,, is the val^e of the exponent, namely r^ - r. Figure 6 shows r^ - r 
as a function of r^ for two representative rays. .. 
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Fig. 5—Dimensionless pressure components for 2 curve of Fig. 
Yield = 10 KT, burst depth = 200 ft. 
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hundreds of feet 

Distance from virtual source, r^ 

Fig. 6—"Phase" of compression wave at time of maximum tension; yield = 10 KT 
burst depth = 200 ft. 
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IV. EXTENT OF CAVITATION 

Having developed the pressure field assuming no cavitation, the next step 
is to determine how this pressure field is altered by cavitation. The basic cri¬ 
teria for the commencement and termination of cavitation were discussed in some de¬ 
tail in the already cited report bearing on the orve-dimensional model of bulk cavi¬ 
tation. The same basic theory is now followed in the development of the extent of 
cavitation in the three-dimensional case. Cavitation will commence at any point 
at which the tension reaches some constant value which might be called the rupture 
or breaking pressure. This pressure is a function pf the fluid being considered. 
There is considerable doubt as to what is the proper value of the breaking pressure 
for seawaterf however, It developed in the one-dimensional case that the exact nu¬ 
merical value of this quantity did not strongly effect the results and so not too 
much attention has been P*id to determining the best value of this constant. In 
the present study two values of this breaking pressure have been used. One is -10 
psi and the other is -100 psi. It is felt that with these two values the range of 
practical realistic values will be represented. 

It might then be thought that the cavitation or rupture process will con¬ 
tinue until the total pressure exceeds this rupture or breaking pressure. However, 
a more careful analysis of the situation indicates that this is not the proper cri¬ 
terion for the termination of cavitation. If we consider that the rarefaction 
front does not have a truly instantaneous rise, theh we can see that as cavitation 
commences the breaking front is eroded away so that the tail of the rarefaction 
wave is lost and that portion of the head above the tension required to maintain 
breaking pressure (or tension) is eroded away and not carried through the cavi- 


3or 
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tated region 4 . Now, cavitation is. clearly a dissipative process and not a regener¬ 


ative one so if at any point, the negative pressure or tension required to maintain 


cavitation must be increased, this increase cannot be physically realized and so 


cavitation must terminate. Now, in the three-dimensional case, not only is the 


pressure at the breaking front being attenuated by the cavitation process, but also, 
due. to the increasing volume of the sphere as the wave moves from the center, there 


is an additional attenuation factor proportional to, the inverse of the distance 


from the source. Hence in the-three dimensional case the; criterion for the termin¬ 


ation of cavitation is that the product of the distance from the virtual source,. 


which is the point of origin of the rarefaction wave,, and the minimum tension or 


negative pressure required to maintain rupture pressure must never increase' in ab¬ 


solute value. Another viewpoint is that the pressure ? r at the rarefaction frosty, 


in the absence of rupture, is described in the acoustic, approximation by 

P r = S/l^ (22) 


where 


S is the (virtual) source.strength; 


r, is distance from the virtual source. 
± 


Now with cavitation it is physically, possible to reduce the effective value of S 


by the rupture process—but-It is impossible to increase .it. Hence, while cavita¬ 


tion is taking place .we have , 


dS !P r ) 

df 1- dt 


^Vincent Cushing; loc. cit. pp. 57-59. 
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When there is no rupture the wave travels aa an acoustic wave, with spherical di¬ 
vergence, and with S equal to a constant, i.e. 

dS jjflV _ n (24) 

dt “ —dv • " u 

After termination of cavitation the rarefaction wave assumes a peak value equiva¬ 
lent to that which it had upon tenaination and thereafter progresses onward with 
the inverse distance attenuation, presumably indefinitely. However, if at any 
time the toted pressure including this attenuation and the reduced rarefaction 
front reaches a value equal to the rupture pressure again., then cavitation will 
again commence and continue until the same criterion as used above for the ter¬ 
mination of cavitation is again reached. 

Now, a priori, it may setm unreasonable tto tftL'k -of two regions of cavi¬ 
tation along a given ray, however, upon finding (the .numerical values for the ex¬ 
tent of cavitation it .was found that in some cases there was a sharp break be¬ 
tween two adjacent points on the -curve representing the termination of cavitation. 
Also, some of the total pressure curves presented in Section III., Figs. 2-4, indi¬ 
cated that the region of cavitation might he rjflpubly-connected. A further investi¬ 
gation of these cases indicated that in the -model used in this study, at least, 
it is certainly possible under some condition to have two cavitating regions a- 
long a .given ray. 

Now to define these criteria mathematically we will pick up equation (20) 
for the total pressure.at the head-of-the rarefaction wave 

r -(r x - r) rt 

P t (z,h) = d I* - M + Gz + P a (20) 



i 
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Nov, when this total pressure is equal to P^, the breaking pressure, which 
it must be remembered is a negative number, cavitation will commence. Thus the cri- 

. .> I* * A , 1. . 1 ‘ . >' . . . ■ . 

teria for the commencement of cavitation for the first time along the ray is 


'-d r i - r) . 11 

r ij 


+ Gz + P - P, = 0 
a b 


(25) 


The negative pressure or tension required to maintain the total pressure at a point 
in the cavitated region at precisely the breaking pressure is given by 


V-V-V-o*-? 


-(r, - r) 


( 26 ) 


The criterion for the termination of cavitation is that the rjP 2 reach a 
minimum value. This must be interpreted as a minimum absolute value since P 2 is 
negative and hence the value of the actual product is a maximum. This condition 
may be expressed by the following equation 


d < r l P 2> 

dr. 


= P 2 - Gr^ cos a 


v-r* 


- 2d cos a 


h 0 (27) 


where a is the angle between the ray being considered and the vertical. Now, if 
the value of r^ which meets this criterion is designated r^ g and the value of P 2 at 
> this value of r^ is denoted by P 2g , then the total pressure in the region below this 
cavitated region can be represented by 

. '-(r, - r) 


de 


p t = “r - + Gz + p a 


h 


( 28 ) 


and cavitation will begin a second time whenever this total pressure becomes equal 
to the rupture pressure. The second-cessation of cavitation will occur when the 
condition given by Eq. (4) is again satisfied. 
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A. Numerical Results 

Using the model described above, a series of numerical results have been 
obtained for various explosions in which both the depth of burst and the explosive 
yield have been varied. In order to obtain these numerical results certain pres¬ 
sures had to be evaluated. These explosive parameters were based on the results of 

previous tests. The following parameters were used} 

, , , . 

P o = 14.7 psi 
G = ,64.4 lbs/cu. ft. 

P fl = 1.38 x loV/ 3 /d psi 
A = 53.5V 2 / 3 ft. 

where V is the yield of. the explosion in kilotons 
d is the depth of the explosion in feet. 

Using these constants the results shown in Figs. 7. through 24 were obtained for 
the extent .of the cavitated region under the conditions noted on each figure. In 
Table 1 the various conditions shown in the figures are listed together with seme 
of the results. 

It can be seen that there are four types of explosions from a cavitation 

% * ■ 

* 

viewpoint. There is one type in which the cavitated region is singly connected 
and above the explosion in the immediate neighborhood of the explosion. Two, 

• there is the type of explosion in which the cavitated region is doubly connected 
and a non-cavitated "ball" exists around the explosion. Three, the case in which 
the cavitated region is singly connected and Is completely below the explosion in 
the immediate region of the.burst point. Lastly, if the explosion occurs at great 
depth, there will be no cavitation at all. A further study was made to determine 
approximately the boundaries between these four regions and it was found that a 
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TABLE I. 


TYPES OF EXPLOSIONS INVESTIGATED 


Case 

Yield 

-KT 

Burst 
Depth 
- ft T 

Breaking 
Pressure 
- psi 

Type of 
Cavitation 

Maximum 
Horizontal 
Extent of 
Cavitation 
- ft. 

Maximum 
Depth of 
Cavitation 
- ft. 

Illustrated 

in 

Figure No. 

1 

10 

200 

-10 

2 

3,360 

460 

7 

2 

10 

200 

-100 

2 

2,400 

445 

8 

3 

1 

1000 

-100 

1 

4,830 

140 

9 

4 

1 

1000 

-10 

1 

7,500 

190 

10 

5 

1 

200 

-100 

1 

2,280 

230 

11 

6 

1 

200 

-10 

1 

3,230 

240 

12 

7 

10 

1000 

-10 

1 

8,210 

310 

13 

8 

10 

1000 

-100 

1 

6,120 

300 

14 

9 

10 

150 

-10 

2 

2,800 

500 

15 

10 

13 

700 

-10 

1 

6,900 

430 

16 

11 

20 

70 

-10 

3 

1,710 

680 

17 

12 

20 

150 

-10 

2 

2,790 

620 

18 

13 

50 

5000 

-10 

1 

17,900 

380 

19 

14 

20 

5000 

-10 

1 

17,100 

300 

20 


Type 1 - Region of cavitation above burst depth in vertical plane of explosion. 

Type 2 - Region of cavitation both above and below burst depth in vertical 
plane of explosion. 

Type 3 - Region of cavitation below burst depth in vertical plane of explosion. 
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Fig. 7a—Extent of cavitated region; yield : 10 KT, burst 
depth s 200 ft., rupture pressure ■ -10 pci. 



Horizontal distance from burst point, ft. 


Fig. 7b--Detail of area around burst point. 


ENGINE 

V* v- 


EKING 


PHYSICS COMPANY 


V 






















































































































































































Horizontal distance, feet 
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Fig. 8—Extent of cavitated region; yield = 10 KT, 

burst depth = 200 ft., rupture pressure = -100 psi. 
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Fig. 9 —Extent of cavitated region; yield = 1 KT, 

burst depth = 1000 ft., rupture pressure = -100 psi 
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Fig. 10—Extent of cavitated region; yield = 1 KT; 

burst depth = 1000 ft., rupture pressure = -10 psi. 
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Fig. 11—Extent of cavitated region; yield = 1 KT, burst depth = 200 ft. 
rupture pressure = -100 psi. 
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Fig. 12—Extent of cavitated region; yield = 1 KT, 

burst depth = 200 ft ., rupture pressure = -10 psi. 













































































































































Horizontal distance, thousands of ft 
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burst depth = 1000 ft., rupture pressure = -100 psi. 
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Fig. 15--Extent of cavitated region; yield = 10 KT, 

burst depth = 150 ft., rupture pressure = -10 psi. 
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burst depth = 700 ft., rupture pressure = -10 psi. 
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Fig. 18—Extent of cavitated region; yield. = 20 KT 
burst depth = 150 ft., rupture pressure = -10 psi. 
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very shallow Explosion products the third type of cavitation, namely, where all 
the cavitation is direotly below the explosion in the region of the explosion. If 
the depth of the explosion is less than the value of X, as defined above,, this type 
of cavitation may be expected. If the depth is greater than X but somewhat less 
than another value dependent to some extent upon t^e yield; then the doubly connect¬ 
ed cavitated region occurs. The largest value of depth for which this condition 
will hold is plotted in Fig. 21 as a function of yield. Fig. 22 shows a plot of 
the dimensionless depth (d/X) which forms this boundary as a function of yield. 

If the explosion takes place at a depth greater than that indicated in Fig. 21, the 
cavitated region will not extend below the burat point'.. If an explosion occurs at 
an extremely great depth,, it is clear heuristically that no cavitation.would exist 
because the negative pressures or tensions , in the rarefaction wave, would never be 
great enough , to produce- a net tension at any point... The-exact buret {teptfc at which 
the existence of a cavitated region would, cease has not. aa. yet been rigorously de¬ 
termined.. However, it appears from preliminary/result a j . that, if P s , .the pressure, 
due to: the compression wave j at the-surface directly above the burst point; is less 
than, about twice the atmospheric pressure;., then there will be no cavitated region. 

B- Discussion 1 .of Results 

It is certainly recognized that, there 1 ! IP no, experimental evidence to sup¬ 
port the finding of cavitation below the burst point".of' the- explosion and, as a 
matter of fact, in reality it is not precisely* expected! that cavitation would exist 
directly 1 below the burst point because the rarefaction.wave;, used in our model, 
would not be transmitted in a linear fashion through the gas: bubble. However, 
there ip a definite possibility that there may be nearly a doubly connected region 
of bhvitation and that these results will furnish an indication of the phenomena 
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W - yield - KT 

-Depth at which cavitation directly beneath burst point ceases as a function 
of yield. 
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which exist under real conditions. In the very shallow explosion It may be expected 
that there would be an annular region of cavitation with a central region in which 
there was no cavitation, then a cavitating ring, and then, of course, no cavitation 
at great distance. In any event, it is felt that the cavitated regions below an 
explosion will dose very qulokly slnoe the bottom surface is rigidly held and will 
not move downward in a way comparable to the upward movement of a slug of water 
when cavitation occurs above an explosion. 


ENGINEERING-PHYSICS 


COMPANY 




43 


V. TRAJECTORY OF WATER HAMMER 

The following mathematical model was developed to obtain an initial 
understanding of the way in which cavitation pulses are generated. For the 

purposes of this first study it is assumed that the layer of water overlying 

1 ., 

the cavitated region is spalled upward with an initial velocity equal to twice 
the vertical component of the instantaneous particle velocity associated with 
the pressure at the front of the sbo :k wave. The time at which spalling takes 

i 

place at any given horizontal range from the explosive source is determined 
from the sound speed of water and the depth of the explosion. It is further 
assumed that this liquid layer is acted upon only by forces due to the accel¬ 
eration of gravity and furthermore that the usual acoustic relationship between 
pressure and partldle-velocity holds. In addition .it ..is assumed that spherical 
divergence, 4/ft, is the only mechanism by which the magnitude of the shock 
wave is reduced with distance from the explosive source. It was found that 

• '» •* V *• 

this very simple- model predicts a cavitation pulse of striking likeness to 
those observed experimentally. 

In setting up, the equations for $be.. above-described model, the time 
required: for the shock to reach a giyen horizontal range, h, when exploded at 
a depth, d, and where the shock propagates, at a velocity, C, is given by 

. t - R/C , (29) 

• 2 2 2 

where* R - h +. d . .The time' required for. the liquid layer to go through its 
trajectory and impact with.,the underlying water is given by 

t = 2u/g , (30) 

Where ii is the initial velocity of the liquid layer; 
g is the acceleration of gravity. 
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Hence the total time required from detonation to water hammer impact is 

' . I . - 

given by 

T « R/C -t 4u y /g , (31) 

where u y is the vertical component of particle velocity associated with the 
Instantaneous pressure at the shock front at a horizontal distance, h, on 
the surface. Using the acoustic approximation, Eq. (31) can be written as 

T - R/C + 4Pd/gpCR , (32) 

where p is the density and P is the peak pressure. We assume that the peak 
pressure at slant range, R, is given by 

P = KW^/R , (33) 

where W is the yield in kilotons; 

K is a constant having a value oif roughly 2 a 10^ when P is in lb/in 2 
and R is in feet. 

Substituting Eq. (33) into Eq. (32) we obtain 

T - R/C + ItfCdW^/gpCR 2 . (34) 

By differentiating Eq. (34) with respeqt to h and taking 1 the reciprocal of 

both sides, we obtain an expression for the horizontal velocity of the travel¬ 
ing hammer, 

V ■ pCgR^/Ih (pfift 3 -8KdW^ 3 )] . (35) 

For increasing h and, hence, increasing R, for fixed d, it is easy to see from 
Eq. (35) that the water hammer velocity approaches the sound speed of water. 
From Eq. (35) we can also determine the horizontal distance at which first 
impact occurs. 
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h = [(eKW^fypg) 2 / 3 -^] 1 / 2 for d s (SKW^/pg) 1 / 2 


(36) 


H = 0 for d > (SKW^/pg) 1 / 2 

Figure 23 shows the horizontal distance to the circle of first impact as a 
function of the depth of burst. From Eq. (36) we see that depending upon the 
depth of blast and size of the bomb the water hammer travels inward toward the 
source of the explosion or outward away from it or simultaneously in both 
directions from the circle of first impact. It is also perhaps interesting 
that there exists a depth of burst for a given yield which maximizes the dis¬ 
tance to the point of first impact. This depth obtained from Eq. (36) is 

d* = (4/3) 3//4 (K/pg) 1/2 W 1/6 , (37) 

and the corresponding maximum first impact point can be determined by substi- 

— # 

tution of Eq. (37) into Eq. (36), to be h = /”2d . A table giving horizontal 
distance of first impact and time of first impact for a few representative cases 


is given below. 

TABLE II. 


Depth 

Equivalent 

Horizontal Range 

Time of First 

Feet 

Weight - KT 

of First Impact - ft. 

Impact - sec. 

70 

20. 

1700 

0.50 

200 

1 . 

1690 

0.51 

200 

10. 

2190 

0.66 

1000 

1 . 

2740 

0.87 

1000 

10. 

3620 

1.12 

5000 

20. 

4800 

2.08 


Using the above equations a number of computer runs were made on an 
IBM 1620 to determine the water hammer trajectory. Figures 24 and 25 show 
typical plots of time of impact and impact velocity versus horizontal distance. 
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Horizontal distance, h 

Fig. 24—Time of impact and impact velocity versus horizontal 
distance for a 10 KT yield at 200 foot depth. 
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thousand feet 
Horliontal distance, h 

fig. 25 —rise of Inpact and lapact velocity versus horizontal 
distance for a 20 K.T. yield at 5000-foot depth. 
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In all of the above cases except the 5000-foot depth shot the closing 
velocity from the circle of first impact is initially unbounded in both direc¬ 
tions. The closing velocity then approaches the sound speed of water in the 
horizontal direction away from the bomb while the closing velocity in the 
direction toward the bomb soon decreases to rather low subsonic values and 
then again rises to an unbounded value as the water hammer approaches a point 
directly under surface zero. In the case of the 5000-foot depth shot the 
inward traveling water hammer did not go subsonic but rather reached a super¬ 
sonic minimum. Other details of closure for the 5000-foot depth shot were 
similar to those of the other shots investigated. 

In <frder to substantiate the use of this simplified model, calculations 

5 

were made which could be compared with experimental data . A satisfactory 
agreement between experimental and theoretical results was obtained. 


5 

E. Swift, Jr., et al, "Underwater Pressures from Underwater Bursts," 
Operation HARDTACK, Proj. 1.1, WT-1606, U. S. Naval Ordnance Laboratory, August 
15, I960. 
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VI. SECONDARY WAVE CONFIGURATION 

Section V has described the trajectory of the point of reloading. This 
traveling impact point is essentially a moving sound source or pressure source, 
and so we nov investigate at least qualitatively the secondary wave patterns 
to be expected from such a moving source. 

Figure 26a schematically indicates the secondary wave pattern soon after 
the .'time of. first .impact. Near the time of first Impact, both the Inward 
traveling: as’ well as the outward traveling water hammers are moving supersoni¬ 
cally (relative to sound speed in seawater). Indeed, analysis with Eq. (35) as 
indicated in Section V indicates that the outward traveling hamper is always 
traveling, supersonically, but at an ever decreasing rate, asymptotically 
approachingrsound speed. Consequently, typically, as Indicated in Fig. 26a 
this supersonically traveling water haunter may be expected to produce a ,r bow 
shock.." Tfiis bow-shook or Mach line is at an angle 0, ds indicated in Fig. 26a, 
described by 

0 ■ arc 'Sin l/M , ( 38 ) 

where Mi is. .the Mach number of the traveling locus of impact. The Inward 
traveling locus of impact also initially moves at a very high Mach ninnber. 

The Mach lines due to these two traveling, impact loci form a secondary shock 
wave as Indicated in Fig. 26a. As long as both traveling impact points are 
moving supersonically, the associated secondary wave will have a shock front. 

In Fig. 26b ve indicate the secondary wave pattern at a somewhat later 
time, in fact at a. time when the velocity of the inward traveling impact plaint 
has fallen well below* sound speed. Most of the secondary wave configuration 
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Locus of 
first Impact 


Secondary- 
shock front 



Fig. 26a— Secondary wave soon after the time of 
first impact. 



Fig. 26b —Secondary pressure wave at a time when 
the inward traveling impact point has 
fallen well below sound speed. 
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still has a shock rise, since these portions of the secondary wave were pro¬ 
duced while both water hammers were still traveling supersonically. But that 
portion of the secondary wave which was generated while the inward traveling 
water hammer, was-moving at less than sound speed, will not have a shock rise 
appearance, but rather will have a gradual rise or precursor type of charac¬ 
teristic. 

In Fig. 26b we indicate that the Mach lines due to the outward traveling 
hammer produce wave fronts which are pointed radially outward. Oh the other 
hand, the inward traveling water hammer—which, we have seen above, produces 
much more Intense secondary waves—produces wave fronts which are traveling 
inward. That portion of the secondary wave configuration which is generated 
by,' the inward traveling waterrhammer tends to drive-or focus secondary hydro- 
dynamic wave energy toward two regions lying • oa the vertical axis passing 
through the explosion point: one point, would appear to He beneath the explo¬ 
sion point; while the other point would appear to be well above the explosion 
point. It is the latter region of concentration which is most meaningful in 
regard to the airblast above an underwater nuclear explosion. 

In regard to the inward traveling water hammer it is interesting to call 
attention to the very strong Mach lir - focusing effects which occur when a 
sound source orrpressure source is accelerating or decelerating through the 
trans-sonic region. This is the "sonic boom" phenomenon which has been widely 
experienced and theoretically substantiated.^ 


R. A. . Struble, C. E. Stewart, E A. Brown nnd A. Ritter, '^Theoretical 
Investigations of Sonic Boom Phenomena, 1 Armour Research Foundation of Illinois 
Institute of Technology for Wright Air Development Center, August 1957* 
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VII. SUMMARY 

The previous one-dimensional analysis of bulk cavitation is being extended 
to the three-dimensional analysis. In the present report, results of a detailed 
analysis of the extent of the region of cavitation and a rudimentary analysis of 
the closure process are presented. 

The detailed analysis shows that the cavitated region is much more exten¬ 
sive than it would appear to be in the traditionally accepted model. In particu¬ 
lar, the cavitated region in our analysis may extend directly under the burst point. 
Although our analysis to date has neglected the gas globe at the burst point (and 
the physical existence of the gas globe will undoubtedly have a significant effect 
on cavitation directly beneath the burst point), it must be remembered that cavi¬ 
tation forms quickly and at an early time when the effect of gas globe is small. 
Thus, while the analytic result must be considered to be approximate directly under 
the burst point, it should be sufficiently accurate in nearby regions. The newly 
discovered shape of the cavitated region may have a significant effect on the water 
hammer since a large proportion of the water hammer energy is concentrated near 
the burst point. Therefore, it is believed that development of the more accurate 
estimate of the extent of the cavitating region, as presented here, will greatly 
improve the overall theory of water hammer. 

In an effort to keep our initial models simple, we have not yet considered 
bottom effects. Consequently models may be considered accurate representations of 
explosions in deep water, (or of explosions on a rigid bottom), but the results 
must be modified to represent effects of explosions in shallow water. The neces¬ 
sary modifications will be developed in the continuing work. 
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The rudimentary analysis of closure presented in this report agrees satis¬ 
factorily with experimental results of large scale explosions and indicates that 
the major effects of the water hammer may be susceptible to approximation with a 
simplified theory for these cases. This preliminary analysis, however, does not 
consider the effect of atmospheric pressure on closure and it is expected that in 
small explosions, this effect will be important and thus require more refined an¬ 
alysis. It is interesting to note that this preliminary analysis assumes that 
the surface will have returned to its initial position at the time of closure. 

Since the theory seemingly agrees with experimental results (at least for larger 
explosions), it appears that distortions of the surface at the time of closure 
must be small. 

The estimate of the secondary wave pattern based on this rudimentary clo¬ 
sure analysis indicates that for most explosions initial impact will occur along 
a ring some (horizontal) distance from the burst point. Closure or water hammer 
will then proceed both radially inward and outward. The outward speed of closure 
will always be supersonic. However, except for very deep explosions, the inward 
speed of the traveling water hammer will initially be supersonic, tnen become 
subsonic and again rise to supersonic as surface zero is approached. As a result, 
the "sonic boom" phenomena will exist with an impulsive release of energy "focused" 
to the vertical line through the burst point. 
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